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In the first order formalism of gravitational theories, the spacetime connection is considered as an 
independent variable to vary together with the metric. However, the metric still generates its Levi- 
Civita connection that turns out to determine the geodesies of matter. Recently, "hybrid" gravity 
theories have been introduced by constructing actions involving both the independent Palatini con- 
nection and the metric Levi-Civita connection. In this study a method is developed to analyse the 

fT^ ' field content of such theories, in particular to determine whether the propagating degrees of freedom 

are ghosts or tachyons. New types of second, fourth and sixth order derivative gravity theories are 

c7^ ' investigated and the so called f{X) theories are singled out as a viable class of "hybrid" extensions 

O^ ' of General Relativity. 
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CN ! I. INTRODUCTION 



Motivated by the oddities discovered by recent astronomical observations, a plethora of modified theories of gravity 

O ' has been advanced during the last years. Their main scope is to reproduce the observed behavior of our universe 

without invoking any undetected entity such as dark matter, dark energy or the infiaton field. Although some of them 

5-H ' may succeed in reaching the goal, some others give rise to unphysical features such as, for example, the appearance 

^ O^ ' of ghosts. 

Among all these modified theories, the most considered are the so-called f{R) theories of gravity, where the action 
y—^ is taken to depend on a general function of the Ricci scalar R\ see [ij for some reviews. Their popularity comes from 
t> I. the fact that they are both theoretically viable and sufficiently simple to study in cosmology and other frameworks. To 
C ■ ■ these ends, the f{R) action is usually varied with respect to the metric tensor in order to produce the field equations 
■^ ,: governing the dynamics of gravitation. This approach is known under the name of metric (sometimes second-order) 
^ir ' variational principle and it is the one first employed by Hilbert to derive the Einstein field equations. However a 
• , second variational method is well known since Einstein formulated it, though for historical reasons [2] it has been 
\l ■ named Palatini (sometimes first-order) variational principle. It consists of an independent variation with respect of 
the metric and the (torsionless) connection. If applied to the Einstein-Hilbert action it gives back the Einstein field 
equations, meaning there are no physical difference for general relativity. However when considered for other theories, 
such as f{R) gravity, it produces completely new dynamical field equations resulting in a different phenomenology; 
see [U, Q for reviews. 

The metric and Palatini approaches have been recently combined^ to give born to a new class of modified gravi- 
;J] , tational theories which has been named hybrid metric-Palatini or f{X) gravity [4|. The action is taken to depend 
linearly on the metric curvature scalar R but nonlinearly on the Palatini curvature scalar R which is modulated by an 
arbitrary function in analogy with Palatini f{R) theories. Some cosmological and astrophysical applications of these 
theories have already been studied and it has been shown that they can both pass the solar system experiments and 
predict a late time accelerated cosmological expansion. For this reason they can constitute a promising theoretical 
explanation for dark energy. The hybrid metric-Palatini approach has then been generalized considering a general 
function of both R and R [Tj ; independently from an earlier work which considered similar ideas for different reasons 
3 . The cosmological evolution of these generalized theories have also been analyzed and late time acceleration can 
be achieved in this context. Moreover it has been proved that this class of theories is dynamical equivalent to a 
nonminimally coupled biscalar theory of gravity and that one of the two scalar disappears in the f{X) subclass. In 
other words, the theory introduces two new scalar degrees of freedom, one in the f{X) case, on top of the metric ones. 
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^ See Refs.l^lgl for other frameworks for unifying the variational principles. 



In this paper we will consider a further generalization of gravitational theories within the hybrid metric-Palatini 
framework. In particular we will take the gravitational action to arbitrarily depend on R, R and Qh — R^'^R^iv 
This generalization is done in analogy with the Ricci square modified gravity which has been largely studied in both 
metric ISIj'lO] and Palatini [11] formulations. Note however that the invariants Qh we are considering can only appear 
within a hybrid metric-Palatini context, and thus were never considered before. Our analysis aims at characterizing 
the flat space propagator in order to verify the presence of ghosts for such theories. It is a well known result that 
Ricci square modifications always leads to a spin-2 ghost called the Weyl ghost [9]. It is thus interesting to perform 
the ghost analysis also for this class of hybrid theories in order to check if also these models are inevitably haunted 
by ghosts. 

The paper is organized as follows. In SeclUwe will introduce the action of the theory and derive the field equations 
discussing how they can be solved in general. In Sec. IIIII we will consider an expansion around Minkowski space and 
compute the field equations in this weak field limit. In Sec. IIVI we will then invert the field equations to calculate 
the flat space propagator. An analysis for several subcases included in the class of theories we consider will be briefly 
carried out in order to check that the results are consistent with previous works, and a health diagnosis will be 
performed for the new theories included in our action. Finally we draw our conclusions in Sec. |Vl In what follows we 
will set K^ = SttG/c* for the sake of simplicity. 

II. ACTION AND FIELD EQUATIONS 

The action we will adopt for our analysis is given by 

S=^Jd^xV^f{R,R,QH), (1) 

where R = g^'^Rfj.n and R = g^'^Rfj.v are the metric and Palatini curvature scalars respectively. The first one is 
composed with the usual Levi-Civita connection F^^ and it is completely determined by the metric 5^1/, 

D ^ j^a pQ I pa pA j^a pA fr^\ 

Ji-ZJi/ — ^ /ii/,Q ^ ^a.v I ^ qA^ ^ii/ ^ fi\^ ay i K'^ ) 

the second one depends on an independent torsionless connection F^^^ forming the following Ricci tensor 

D ^^ pCK pQ I pa pA f^OL pA fo\ 

J^fif — i ;^i/,Q ^ f^a.f ' ^ qA^ fiv ^ ^A^ ail ■ K"^ ) 

The function / in action ([T]) arbitrarily depends on R, R and 

Qh = R^'R^.u , (4) 

which will be denoted as the hybrid Ricci square invariant. Note that, though i?^!/ is in general asymmetric, only its 
symmetric part enters the action ([1]). 

Action ID) generalizes the theory considered in 7] where the hybrid Ricci square invariant was absent. As shown 
in [T], |8| such a theory is dynamically equivalent to a nonminimally coupled biscalar theory. Here we undertake the 
study of viability of such a biscalar theory. It is interesting to consider also theories with other higher order invariants 
in the hybrid metric-Palatini context. Although several new invariants can be built within this framework , we decide 
to work only with Qh since it is the simpler among the higher curvature invariants not appearing in pure metric and 
Palatini gravities and thus never considered in the previous literature. 

Variation of ([T]) with respect to the metric produces the following field equations 

+ 2fQR'^K, + in, (/gi?^.) + ig^.V^V^ (/q^"^) - VaV(. {f.^R^)) = ^^T,, . (5) 

where /.a, /a and /a are the derivatives of / with respect to i?, R and Qh respectively. Round brackets between 
two indices denotes symmetrization and D, — g'^'^V^Vj/. The energy-momentum tensor on the right hand side of 
([5]) has been derived adding a standard matter action to ([T]). Variation of action ^ with respect to the independent 
connection F^^, yields, after some manipulations, the condition 

Va (V^/^ff^'' + V~9f,QR^'') = , (6) 



where V^ denotes the covariant derivative of f ^y. In order to solve this equation a new metric g^^ imphcitly defined 

by 

V^gr = V~9 [fM^" + f^QR"") , (7) 

emerges in Palatini theories. If the right hand side of (O is independent of the connection F^^ a solution is given by 

r^. = 1?^" id^9a. + d,g,^ ~ d^g^,) , (8) 

or, in other words, the independent connection is given by the Levi-Civita connection in terms of the metric 5^^. 
However in ([7]) the two functions fh and /a still implicitly depends on F^^. Within the purely Palatini framework 

one usually employs field equations ([5]) in order to algebraically relate R and Q to the matter fields [lli| . In this way one 
can substitute these invariants with matter fields wherever they appear in the field equations. In particular replacing 
them inside the condition ([7]) allows to solve for the independent connection as in (|S]). In our case however this cannot 
be done in general since in the field equations ([5|) derivatives of Qh explicitly appear and thus an algebraic relation 
with matter fields is impossible to achieve. Although this prevents finding general solutions of the field equations in 
a simple algebraic way, we will see in the next section that in the weak field limit needed for the ghosts analysis the 
right hand side of (O will be independent of F^^ at first order in the perturbations. 

The explicit form of the (inverse) metric g'"^ can be read off from Eq.© by taking its determinant. Defining 



rA*" — f . ^M!^ 



1,0""^ + f,QR^'' ' (9) 



and 1/r = det(r^'^) we find 



g"" = ^^ r'^" . (10) 



III. WEAK FIELD LIMIT 

We now consider an expansion around Minkowski spacetime. The metric tensor will be then given by 

.9^11/ = r/i-Li' + hi_,y , with |/ip^| < 1 , (11) 

where 77^1, is the Minkowski metric. All the standard properties of the weak field limit approximations will hold, in 
particular h^^i, is Lorentz covariant and can be treated as a Lorentz tensor on Minkowski spacetime. Up to second 
order in h^,, we have g^" = -q^" — h^'^ + o{h?) and 

g = -l + h + o{h'^), (12) 

where we define h = rj'^'^hfj.^, and denote with o{h'^) second and higher orders in h^^. 

At this point we need to expand / and its derivatives. Since these are functions of R, R and Qh we need first to 
know which order in /i^i, are these invariants. The metric curvature scalar R is of course of order one and has no 
zeroth order term as it happens in GR. The Ricci tensor formed with the Palatini connection do not depend explicitly 



on the metric. However if one had a solution of f ^^^ in terms of 5^1/ this problem would not be present. In what 



follows we will assume that also i?'"^ is of order one in /i^ii^ and that its zeroth order term vanishes. This ansatz will 
be verified at the end of the calculations where we will indeed show that this is true. With this assumption we have 
R ^ o{h) and Qh ^ o{h?). The expansion of / (and thus of all its derivatives) will be given by 

/(i?, R, Qh) = ./(°) + f^'j^R^'^ + f%^R^'^ + o{e) , (13) 

where /^°-', /^ and / - are the constant values of the functions evaluated a.t R — R — Qh — and R'^^\ RS^^ are 

the first order terms of R and R respectively. 

In what comes next we will make another ansatz, namely / ^ a = or, in other words, the second derivative of / 

with respect to R has to vanish when evaluated dX R = R = Qh = 0. This assumption is needed in order to solve 
^ in the weak field limit approximation, i.e. it allows r**" to be F^^, independent in this limit. Note that this puts a 



constraint on the general expression / can assume. If / is independent or linear in R this assumption is automatically 
satisfied. However this is true for more general functions / and thus we will still consider a general R dependency. 
Within these assumptions the expansion of r^"^ is given by 

rt^'^ = f<^l^r]t"' - f'^l^ht"' + /^°' i^fi^^" + /(°)i?(i)^'' + o(h^) . (14) 

One can immediately notice that now r^'^ is F^^ indepedent up to first order in hf^^. Thanks to ([H)l we can also find 
the expansion of 1/r which is given by 

/ 4J-(0) , .(0) 

- = f'f \-^-h + -^^^(Sr^^W 1 + oih^) . (15) 

We are now ready to expand g^'^ given by pUj) . Up to first order in h^i, we find 
where 



r^^ff; iv'"'' - h^^" ) + om , (16) 



qAO) ^ ^(0) „(0) 

h'"' = h'"' + — 5^L_Ai?(i)^A''' + ^i?(i)A''> . (17) 

^•^,A ■^R 

The inverse of 5^" is given by 

5m- = /°^"' {v^.. + K-^) + o{h^) , (18) 

and one can indeed verify that g^°'gau — 5'i + o{h^). We can now give the solution of F^^ up to first order in h^i, 
using ([5]). The result is 

r^. = \v^" [dX- + 9. V - 5- V) + o{h^) ■ (19) 

Note that the constant factor in front of P^ is completely unimportant since cancels out in Ff^^ and then will not 
appear in the field equations. This is in agreement with the fact that g^u can always be rescaled by a constant factor 
without altering the physics. 

The expression (|19p for F^^, allows us to write down the expansion of R^^ as 

hi^ = Ri^J + o{h^) = a.a(>^) - ^d^d^h - ^Dh^, + o{h^) , (20) 

where of course h = i^^^h^y. One can now realize that our previous assumptions that R^^ has a vanishing zeroth 
order term is indeed verified. We can now expand -R/iJ in terms of h^^. Defining 

A = 7^(0) '^ - --d 5 = -^, (21) 

we find 



We have now all the ingredients to expand the field equations ([5]) up to first order in h^i, . Following [1^ |13| we can 
write them in the form 

^ {Vf.. + Zip.) + a{U)Uh^, + 2b{U)d^d(^hl) + c(n) {rj^.d^dph''^ + d^^d^h) 

+ d{U)r^^,Uh + £Ma^a^a^a^/j"^ ^ -2^2 y^_^ ^ (23) 



where a(n), 6(n), c(n), d{0) and e(n) are functions of □ = rj^^^dp^d^- In these equations Z*^"^ represents the cosmo- 
logical constant and will not influence the upcoming analysis. In our case the functions appearing in (|23p are given 

by 



a(n) = 


f(o) , .(0) f(0)^n2 

-IR +•^,7? "^Q T ' 














(24) 


b{u)-. 


=-/!];,^-/5,^+/°'f°^ 


(25) 


c{a) = 


^/;?'+/l^-2(/,a+4/i+/g))n+ 


C(6A + B) + /2)(2A+f); 


tf, 


(26) 


d{a). 


-Z]?^ -/°^ + 2 (/(];), + 4/°)^ + /g))n - 


/l(6A + i^) + /S^(2A+f); 


^^ 


(27) 


e(n) = 


^2(^S,u<,+/S0°~ 


/S(6A + B) + /2)(2A+f)' 


n2. 






(28) 



We can immediately notice that the following relations hold 



0. 



(29) 



At first order these constraints must be true in any gravitational theory preserving difFeomorphism invariance inasmuch 
as they directly follow from the Bianchi identities p^.ll3j. Note also that in the GR limit f -J — >■ 1, / - -^ 0, / - — >■ 



(0) 



(or equivalently f^j^' -^ 0, f 



f(o) 



1 f'"' 



0) we recover the right values a 



1 and e 



0. 



With the field equations expanded at first order in /i^^ we can now compute the flat space propagator and check if 
ghosts appear in the theory. 



IV. FLAT SPACE PROPAGATORS 



In this section, we will analyze the field content and viability of the theories contained in our general class of 
actions. We adopt the formalism presented in Ref. ^13] . that was originally introduced by van Niewenhuizen in Ref. 
|14| . employed in Ref. [15] to study the class of /(i?) models and recently generalized in Ref. [I4I for general metric 
theories of gravitation involving arbitrary combinations of curvature invariants at an arbitrary order in derivatives. 
The main result is that the propagator 11 for metric fluctuations around flat background can be written, in the Fourier 
space, as 

''^^ a{-k^) ~ a{-k2) ~ 3c{-k^) ' ^^°^ 

where V^ picks up the spin-2 and V'^ the scalar modes of the fluctuations. 

Let us start by considering the class of f{R, R) theories. Without significant loss of generality, we may assume that 
f j^ + / - = 1, so that we recover Einstein's General relativity for the graviton propagator at the infrared limit. The 

sum may converge to another positive constant than unity in this limit, but taking that into account amounts to an 
overall rescaling which is irrelevant for our purposes. We readily see that then 



a — —b — 1 , c — —d — 1 
The propagator we then obtain is 
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where the GR propagator is given by 
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n 



vZk- 



GR 



1 



V' 



I2f%k- 



zV 



^Af%k- 



(31) 



(32) 



(33) 



A subtlety here is that though the scalar propagator comes with the wrong sign, it does not imply a ghost but rather 
cancels the unphysical longitudinal degree of freedom contained in the spin-2 propagator V^ [Ij, (Tj] . Let us first 
check that the result agrees with the previously known results in other appropriate limits and then proceed to study 
the new generalizations. 



A. Metric f{R) models 



f(0) 



In the pure metric f{R) case, / - — A = and we have [15| 

n/(fl) = Hgr + —. ^-— ^T'" . (34) 

Thus we have an extra scalar degree of freedom, as we expect since the f{R) models are known to be equivalent to 
Brans-Dicke theories with a vanishing parameter ujbd = 0. The mass of the "scalaron" is iri^ = (3/)^^"^, and as 
long as f"{R) > the theory is stable, otherwise a tachyonic mass spoils the stability around Minkowski space. This 
agrees with the stability condition found in [lfl| and generalized in Q- 

B. Palatini f{R) models 

It is well known that the Palatini- type f{R) models are equivalent to Brans-Dicke theories with the parameter 
ujbd = —3/2. This particular value corresponds to vanishing kinetic term of the field, which thus is nondynamical. 
Therefore we expect that no additional scalar degree of freedom should appear. As previously, we may assume that 
f ■ = 1, and we have now of course that f pV, = P K = f\' = 0. Hence, 

^f(R) = ^GR , (35) 

confirming our expectation. 

A comment here has to be made for the sake of generality. Recall that we assumed the condition /ki = which 
actually restrict our analysis in the nonlinear Palatini case. However it easy to realize that Palatini f{R) gravity 
adds no degrees of freedom from some considerations in absence of matter fields. In vacuum the f{R) gravity field 
equations are 



1 
2 



f^iiR^iu - 77.9^1./ = , (36) 



{V^f,R9''')=0- (37) 



Once the function / is prescribed, the trace of (|36|) will give an algebraic equations for R whose solutions will generally 
he R = Rq = const. This implies that /^ will be a constant too and then equation (|57)) will reduce the independent 
connection to Levi-Civita. The field equations ([36|) will then become the Einstein field equations with a cosmological 
constant in vacuum which is a well known result of Palatini f{R) gravity. The propagator of the theory will thus be 

equal to the GR one, as shown in p5l) . even if the /ii = constraint is relaxed. 

RR 

C. Hybrid f{X) models 

The lagrangian of the form R + f{R) results in field equations where the correction to Einstein gravity are controlled 
by the failure of the GR trace equation X = R — k^T. They belong to the algebraic class of scalar-tensor theories 
which single out the unique interpolation of the metric and Palatini f{R) theories by requiring that the scalar field 
has an algebraic expression in terms of the trace and the Ricci scalar, (j> — (j)(T,R) [ITj- It was mentioned in Ref. 
[a| that this produces the peculiar f{X) type of field equations, and their implications and phenomenology has been 
recently studied in contexts of Solar system constraints, cosmology and wormholes [J]. It was already remarked in 
[17| that in Ricci-flat spacetimes the f{X) theories share the properties of Palatini f{R) theories, which in vacuum 
reduce to GR with a possible cosmological constant. Therefore it is not a surprise that we find no new propagating 
degrees of freedom in Minkowski vacuum, 

H/(x) = Hcfl . (38) 

Interestingly though, this class of theories is not equivalent to either of the previous two cases, since when one considers 
curved spacetimes, a new scalar degree of freedom appears. In this sense, the f{X) gravity is a quite minimalistic 



Again we notice that the constraint /ki = prevents a complete analysis also in this case. However we can repeat 

the 

^A"' - 7i9^iyR + f,RRf^<y - T^Qt^yJ = , (39) 



scalar-tensor extension of GR, as the scalar propagates only in the presence of curvature. For the viability criteria for 
these models in regards of cosmology and Solar system, we refer the reader to Q . 

f(o) 

RR 

the same argument above and consider the f{X) field equations in vacuum 

1 „ . - 1 

2 

/^/^5'^"^)=0, (40) 

where now the Lagrangian is i? + ,f{R) and the function / only depends on the Palatini curvature scalar R. If we 
now take the trace of (j39| we obtain, once the function / is prescribed, an algebraic equation relating R to it!. The 
function /^ will generally become a function of R and equation (|40l) will reduce to its (metric) f{R) correspondent. 
The propagator analysis will then equal the metric f{R) gravity one meaning that in the end we will find a scalar 
propagating degree of freedom in addition to the graviton. As noted in [j| the theory is thus quite different from 
Palatini f{R) gravity presenting a richer phenomenology. Note that when the metric Ricci scalar vanishes, as in 
Minkowski spacetime, the trace of equation (p9|) will again give the solution R — Rq and the scalar degree of freedom 
will disappear in agreement with ()38p . As said above the f{X) gravity phenomenology differs from the f{R) gravity 
only in curved spacetimes. 

D. The hybrid f{R,R) models 

The generalized hybrid Ricci scalar theories were introduced in [3, Q and found to have qualitatively different 
properties compared to the more restricted class of f{X) models described above. In particular, the f{R,R) were 
shown to be equivalent to a class of biscalar-tensor theories. It is indeed seen from our formula (j32|) that these theories 
have an extra V^ spin-0 propagator with a double pole, corresponding to two propagating scalar degrees of freedom. 
From the formula, we can easily deduce the masses of these scalar fields. We get 



(0) 



,R 



18 / .(0) 



(/ 



RR 



(/a + <^±^), (41) 



where we have defined for convenience 



\ 



(/. 



RR^^^,rr) fio) ■ ^ > 

,R 



We note that the scalar particle with mass squared m'^_ corresponds to the scalaron appearing in (j34[) in the limit 
of pure f{R) gravity, but in general now has a shifted mass. The other scalar is a new particle that occurs due to 
nontrivial dependence upon R, and unlike in the case of f{X) gravity, it propagates also in Ricci-flat spaces. The 
condition that neither of the scalars has a tachyonic instability, is given by 

/(|>0, and /W +4/(;;)^-5>0. (43) 

The residues at the two poles corresponding to these masses are 

r± = ^-^ ^^ . (44) 

In order for neither of these scalars to be a ghost, we should have both r_|_ > and r_ > 0. The second condition 
would require that 

f.RR + ^f%-S<Q, (45) 

in contradiction with (j43l) . It seems then that we cannot avoid both tachyons and ghosts in this theory. 
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f(0) 



Again we recall that we have assumed /aa = in our analysis. It might be that relaxing this constraint allows for 

viable f{R,R) theories. This would anyway going to depend on the matter degrees of freedom since R can only be 
replaced in the field equations solving the trace of ([5]). Besides the problems we already mentioned which disappear 
with the /a A = ansatz, a complete analysis for this class of theories cannot be performed without explicitly specifying 
the function / and a matter action. 

E. The hybrid Ricci-squared f{R, Q) theories 

Let us finally consider the Qij-invariant. For simplicity, we restrict to models here without nonlinear dependence 
on the metric Ricci scalar; it is easy to see that this does not affect our conclusions essentially. Basically the graviton 
propagator acquires its structure from the function a(n) in (1241) . and now only the higher-derivative term Qh modifies 
it. We can arrange the result for the propagator in the form 

Han 3/g)(l + |/°'P) 

n ~ -■, = \ '^ ^ '^ ^ p° (46) 

The sixth order theory we have at hand has a modulated graviton propagator which adds two extra poles. In addition, 
there appears a scalar propagator that has five poles. This is in a quite drastic contrast with respect to the metric 
Q-theory which contains only one additional spin-2 particle and features fourth order field equations. We need not 
analyze in detail the properties of the new degrees of freedom here, since it is obvious the theory as such is seriously 
haunted by ghosts and thus not physical. 

V. CONCLUSIONS 

We considered theories of gravity in the generalized hybrid framework where, besides the independent Palatini 
connection, the metric Levi-Civita connection is also allowed to enter into the fundamental action. In particular, we 
applied the formalism of |12h15J to uncover the classical and quantum stability of the theories, i.e. to find out which 
of the theories may be free of ghosts and tachyonic instabilities. 

To this end, we started with an action involving three types of curvature invariants constructed from the metric and 
the independent connection, and brought the general theory into a form of an effectively purely metric theory at the 
relevant limit. The generalisation of the analysis to other types of invariants is not entirely straightforward because 
an algebraic solution for the independent connection in terms of the metric and the curvature invariants may not be 
generically available, but instead one would need to supplement the system with differential equations of motion for 
the connection. 

For the simplest theories we reproduced the well known results: the metric f{R) is equivalent to a scalar-tensor 
theory and stable given J"{R) > 0, whilst the Palatini-type f{R) theory carries no extra degrees of freedom. Theories 
involving both Ricci invariants show quite interesting properties: the f{X) theories with only linear dependence upon 
the metric R, are viable and in fact reduce to GR in flat spacetime, whereas more general theories with nontrivial 
dependence upon both of the Ricci invariants inevitably contain instabilities. They introduce two new scalar degrees 
of freedom, where either one of them is a ghost or the other has a negative mass-squared associated to its excitations 
around a flat background. 

It is well understood that in the metric formalism, the f{R) theories are a special class of viable GR extensions 
[18| - what we have found here is that the f{X) gravity seems to have a similar unique place amongst the space of 
all theories within the hybrid formalism. Though we have, for the technical reason stated above, explicitly derived 
the propagator only for theories involving the Q besides the Ricci scalar invariants, from the general structure of 
the equations we expect that similar conclusions would hold quite generically for all the hybrid theories. Discovering 
either exceptions (as the Gauss-Bonnet term in the metric formulation) or a robust no-go theorem, is a challenge for 
future studies. 
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